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Identification of Signal Bias in the Variable Flip Angle 
Method by Linear Display of the Algebraic Ernst 
Equation 

Gunther Helms, 1 * Henning Dathe, 2 Nikolaus Weiskopf, 3 and Peter Dechent 1 



A novel linear parameterization for the variable flip angle 
method for longitudinal relaxation time 7^ quantification from 
spoiled steady state MRI is derived from the half angle tangent 
transform, t, of the flip angle. Plotting the signal S at coordi- 
nates x = St and y = Sh, respectively, establishes a line that 
renders signal amplitude and relaxation term separately as y- 
intercept and slope. This representation allows for estimation of 
the respective parameter from the experimental data. A com- 
prehensive analysis of noise propagation is performed. Numeri- 
cal results for efficient optimization of longitudinal relaxation 
time and proton density mapping experiments are derived. 
Appropriate scaling allows for a linear presentation of data that 
are acquired at different short pulse repetition times, TR « T1 
thus increasing flexibility in the data acquisition by removing the 
limitation of a single pulse repetition time. Signal bias, like due 
to slice-selective excitation or imperfect spoiling, can be readily 
identified by systematic deviations from the linear plot. The 
method is illustrated and validated by 3T experiments on phan- 
toms and human brain. Magn Reson Med 66:669-677, 2011. 
©2011 Wiley-Liss, Inc. 
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The degree of longitudinal relaxation time (TJ weighting 
in spoiled gradient echo MRI (1) is imposed by the par- 
tially saturated steady state that is determined by the 
interaction of flip angle and T x relaxation during the rep- 
etition time, TR, (2,3). This leads to the well-known tran- 
sition from proton density weighting to T a weighting. 
The flip angle dependence is exploited in the most com- 
mon quantitative MRI application that estimates TV by 
variable flip angles, VFA, (4-8) or its dual flip angle ver- 
sion. After scaling by the sine and the tangent of the flip 
angle a linear relationship between the corresponding 
signals is established (5) allowing for computationally ef- 
ficient linear programing. Recent research on the VFA 
method focused on noise optimization (9-12) and espe- 
cially bias reduction (13-15). 

We have recently proposed an exact algebraic formula- 
tion of the Ernst equation for the spoiled steady state sig- 
nal (16). A simple rational function is obtained after suit- 
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able nonlinear transformations in flip angle and TR/T a . 
For small flip angles or short TR, these transformations 
conform to proportionality with a third order error. Thus, 
the algebraic signal equation is particularly useful for fast 
low-angle shot (FLASH) applications (17). Here, we show 
that the algebraic Ernst equation provides a highly con- 
spicuous linear display of VFA data. Intercept and slope 
represent signal amplitude and relaxation term, respec- 
tively. Thus, deviating signals, due to residual transverse 
coherences, can be easily identified by visual inspection. 
This article gives a full account of the underlying theory 
and experimental validation by measurements at 3 Tesla. 

Theory 

The common Ernst equation applies to the situation 
where the transverse coherences are irreversibly 
dephased, that is, perfect "spoiling." For given flip angle, 
a, and TR, the signal is given by 

c , • 1 - eorp(-fliTR) , . 1 - E, 
£> = A sma ; — - „„, = A sin a - 



1 — cos a exp(— f^TR) 



1 — cosa E-i 



[1] 



Here and in the following, we use the relaxation rate, 
R-i = l/TV, and E x = exp{— J? a TR) to simplify the notation. 
A = t4(TE) denotes the amplitude of the gradient echo at 
the echo time, TE, after excitation by a = tt/2 under fully 
relaxed conditions (TR ;» T a ). As first suggested in (5), 
the Ernst Eq. 1 can be rearranged into 



S(a)/sin <x = A[l- E x ] + E t ■ S(a)/tan a, 



[2] 



to derive i? a and A from the linear relationship between 
S/sin a and S/tan a. 

We recently employed the half-angle tangent substitu- 
tion of the trigonometric functions for arbitrary flip 
angles below tt (16) 



t = 2 tan a/2 



[3] 



and an analogous hyperbolic tangent substitution of the 
relaxation term 



Pi = = 2fGmi](fl 1 TR/2) 

1+ Hi 



[4] 



to transform the Ernst Eq. 1 into a low order rational 
function in t and p a : 



S(t) 



At 
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© 2011 Wiley-Liss, Inc. 



669 



670 



Helms et al. 




FIG. 1. The algebraic signal equation and its linear plot, a: The normalized signal plotted for three values of exp(-R|TR): 0.9 (short 
dashes, R.,TR = 0.10, = 0.105); 0.6 (long dashes, R,TR = 0.511, P1 = 0.5); and 0.2 (solid, fi-|TR = 1.61, p A = 1.33). Note that a = it 
is projected to infinity, b: The same examples as in a are shown after linear parameterization. The y-intercept corresponds to S/t = A 
toward a = 0. The x-intercepts correspond to the data toward a = it. Ernst conditions are met at SWA = 0.5 (dashed horizontal line). 
The line through the origin represents the signals measured at t = 1. As seen in a, these fall beyond the Ernst angle (short dashes), right 
onto it (long dashes), and below it (solid line). 



The symbol p 1 must not be confused with the proton 
density, which is inherent to the signal amplitude. This 
function is shown in Fig. la. By the transformations [3] 
and [4], the range of sequence parameters is transformed 
from (0 < a < tt) and (0 < R t TR < oo) to (0 < t < cm) 
and (0 < p! < 2). Toward zero, the transformations 
[3] and [4] conform to a and i?iTR, respectively, with a 
third order error. The Ernst angle in terms of t 
is determined by the root of the derivative with respect 
to t: 



t e = A/^Pi = 

This yields the maximum signal 
S E = Ate/2. 



[6] 



[7] 



Since Eq. 5 can be normalized with respect to S E and 
t e , it follows that 



S E Vte t 



showing that there is a linear relationship between S(t)t 
and S(t)/t. After division by t and multiplication of the 
denominator, Eq. 5 can be rearranged as 

S/t = A- S ;t/2 Pi . [9] 

The convenient features of this linear plot of the 
ordinate y = S(t)/t against the abscissa x = S(t)t are 
illustrated in Fig. lb: The negative slope depends only 
on the relaxation term p a . The intercept of the ordinate 
is identical to the signal amplitude, being the limit of 
S(t)/t for t^0. The intercept of the abscissa corre- 
sponds to an infinite t or a = it, The maximum signal 
at the Ernst angle is always found halfway between the 
intercepts since Eq. 7 can be rewritten as y = S e /t e = 
A/2. 

Any fractional signal level p, given by pA = S(t)/t = 
pA, corresponds to a multiple of t e : 



-2pj = T E 



[101 



If N measurements are to be performed to yield AT 
equidistantly spaced p[i), excluding the values 0 and 1, 
p[i] are given by i/[N + 1). For these, Eq. 10 yields 



N + l-i 



-Te, 



fill 



which in turn can be converted into the corresponding 
flip angles. Compromises can be imposed on equidis- 
tantly sampling the actual signal behavior by N measure- 
ments, if the MR system restricts the choice of flip 
angles to integer values. Note that because t/te = 
(tn+i-j/te)" 1 , the corresponding signals should be equal. 

For N signals S t at arbitrary flip angles ay, the relaxa- 
tion parameter p a is calculated from the regression for- 
mula of the slope of y over x. 



1 (PS™) -NEyffi? 
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with Tj obtained from a,- by Eq. 3. This is then trans- 
formed into 



Ti = — TR/ log Ei = TR / log 



2 + Pi 
2 - Pi 



[13] 



The amplitude is given by the standard formula of the 
y-intercept: 

The case N = 2 yields the dual flip angle solutions 
presented in (16). 

Propagation of Signal Noise 

For any given set of independently measured signals S, 
the propagation of signal noise (a s ) into the p t maps can 
be directly from calculated Eq. 12 as 
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Table 1 



Numerically Derived Settings for Minimal Noise Propagation 


N 


a 2 p1 x A%1 2 / aV 


T-IopAe 


T 2opt/ T E 


<t 2 a x P1 2 / aV 


TIAopAe 


T2Aopt/l"E 


2 


4.000 


0.4142 


2.4142 


5.6133 


0.4903 


3.1461 


3 


2.9313 


0.3718 
0.3718 


2.2108 


2.9886 


0.4633 
0.4633 


2.7927 


3 


2.9313 


0.4520 


2.6939 
2.6931 








4 


2.0000 


0.4142 


2.4141 


2.8067 


0.4904 


3.1721 






0.4142 


2.4141 




0.4903 


3.1495 


4 








2.0841 


0.4456 
0.4456 
0.4456 


2.6310 


5 


1.6532 


0.3898 
0.3898 
0.3898 


2.2885 
2.2892 


1.6203 


0.4323 
0.4323 
0.4323 
0.4323 


2.5303 


5 


1 .6532 


0.4370 
0.4371 


2.5664 
2.5656 
2.5633 








6 


1.3333 


0.4142 


2.4142 


1.8711 


0.4903 


3.1483 






0.4142 


2.4142 




0.4905 


3.1471 






0.4142 


2.4142 




0.4904 


3.1471 


6 








1 .3360 


0.4217 

0.4217 

0.4217 

0.42165 

0.42165 


2.4590 



formalized unit-free variances and normalized t as defined in (16).The left side of the look-up table pertains to noise propagation into 
the normalized maps of pi (relaxation), the right side to normalized maps of A (amplitude). 
Variation in the last digits is due to the cut-off condition of the algorithm and rounding. 

Note that the unit-free normalized variances for the symmetric schemes (boldface) are proportional 1/A/, because the dual angle settings 
are replicated. 



and correspondingly for the A maps from Eq. 14. This 
determines the signal-to-noise ratio (SNR) of the parame- 
ter maps. For N > 3, the polynomial order is too high to 
be solved analytically. 

Since the algebraic formulation is computationally effi- 
cient, the optimal settings for t/t e and the corresponding 
minima in the normalized variances (Table 1) were 
obtained by gradient search for a grid of starting values 
programed in Mathematica 4.0 (Wolfram Research 
Europe, Long Hanborough, Oxon, UK). To estimate p a an 
even N required a repetition of the optimized dual angle 
experiment with the variance proportional to 1/N. 
Schemes with uneven N were less signal-to-noise effi- 
cient, but also required only two flip angles, below and 
above the Ernst angle. For A, a similar dichotomy was 
observed; however, the errors decreased with the number 
of points measured below the Ernst angle. Thus, N - 1 
low flip angles yielded the absolute noise minimum, 
which decreased stronger than 1/N. 

If TV is to be mapped with a certain number of meas- 
urements for a given TR and a target value of T lt p a and 
t e should be calculated from Eqs. 4 and 6, so the recom- 
mended values of the single ti can be obtained from the 
table and the respective values of the flip angles are 
obtained from Eq. 3. 



Noise Propagation with Correlated Errors 

The image noise is scaled individually along abscissa 
(<tst) and ordinate (o- s /t). The errors are equal for t = 1 
or a = 53°. At smaller flip angle (t — > 0), the y-error is 
enhanced and the x-error is diminished, and vice versa 
for t — > oo {tx — > it). 

In a function F of two variables, x and y, the statistical 
errors (given by the variances var x and var y) are propa- 
gated by the respective partial derivatives dF/dx. If the 
errors in x and y are not independent (i.e., correlated), 
their covariance has to be taken into account by using 
the general law of Gaussian error propagation. 

varF(x,y) 

i auia o!7/o \( varx cov(x,y)\(dF/dx\ 

= {dF/dx dFjdy)[ U^/t ' I 16 ' 

\cov(y,x) vary J\oF/ayJ 

Here, F is the residual F[x,y) = y — [A — x/2p a ) with 
x = S x t and y = S/t. Thus, the covariance is identical 
to the image noise variance o- s 2 and independent of t. 

dx dv 

cov(x,y) =— as^as = a|. [17] 

This is in stark contrast to the common linearization 
of Eq. 2, where the errors and their covariance increase 
strongly for small flip angles by 1/sin a or 1/tan a (10). 
Equation 16 yields 
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var(y -A + x/2/pl) = (~ + ^) <j% [18] 

The expression in brackets denotes a T-dependent scal- 
ing factor of the image noise. It can be cancelled by 
imposing suitable weights w[t) onto the square residues 

in the linear least squares objective function. These 
weights are maximal at t e and decrease towards the 
intercepts with the x- and y- axes. 

Short TR 

If i? a TR <S 1 then p a can be approximated and Eq. 9 
arranged as 

S/t = A-(S' t/2TR) ■ T t [20] 

Under this condition, St/2TR can be used as abscissa 
to incorporate varying values of TR into the regression. 
Estimates of T a are obtained directly as the (negative) 
slope 

T = 2 W£, Sj/TR, - J2j WTR,- Sthi 
(EiSm/TR^-NEiSfTl/TR? 

This is the generalization of the formula for the dual 
angle experiment (18). 

Flip Angle Inhomogeneities 

If the size of the imaged object is of the order of the ra- 
dio-frequency (RF) wavelength, the transmitted RF field 
and thus the flip angle is not homogeneous across the 
object. The (local) signal in 3D MRI is then determined 
by the spatial distribution of the actual flip angle 

a = fix) a nom [22] 

The factor fix) denotes the bias field describing misca- 
libration, flip angle inhomogeneity, or the profile of slab- 
selective excitation. It has to be determined by suitable 
techniques (19) and taken into account when calculating 
Tj. Otherwise, it will give rise to systematic errors in p a 
and A. In good approximation, t can be replaced by a (in 
radian) in Eqs. 9 and 12 for small flip angles as usually 
used in FLASH (16,17). At 40° flip angle, t is underesti- 
mated by 4%. Apparent estimates of p 1 and the corre- 
sponding relaxation time. 

T lapp w TR/ Plapp « ^/(x) 2 . [23] 

are then calculated from the nominal flip angles in Eqs. 
12 or 21. It can be corrected post hoc by division of the 
flip angle bias field fix) 2 (17) or by applying a unified 
segmentation algorithm (20) for bias correction (21). 

In slice-selective 2D excitation, the signal represents 
the distribution of flip angles across the slice profile, 
which is difficult to describe analytically. This also 



applies to the residual unspoiled transverse magnetiza- 
tion. For such problems, appropriate phenomenological 
correction methods of have been suggested (14,22). 

MATERIALS AND METHODS 

Experiments were performed on a 3T clinical MR system 
(Magnetom Tim Trio, Siemens Healthcare, Erlangen, Ger- 
many). A transmit-receive knee coil was used for phantom 
measurements, and an eight-channel receive-only head coil 
(In vivo, Gainesville, FL) with body coil transmission for 
human studies. Three-dimensional FLASH MRI with non- 
selective excitation was performed to avoid additional bias 
related to the slice profile (see below). This customized 
sequence (21) allowed control of the phase increment used 
for RF spoiling. One healthy adult subject was examined af- 
ter giving informed consent as approved by the ethics com- 
mittee of the Gottingen University Medical Center. 

Experiment 1 : Nonselective Excitation 

To illustrate nearly ideal behavior, experiments were 
performed on 6% agar (by weight, Sigma-Aldrich Sverige 
AB, Stockholm, Sweden) in a spherical phantom of 10 
cm diameter. T a = of 2.41 ± 0.36 sec and T 2 = 20 ± 1 
ms were determined on two different occasions by sin- 
gle-slice 5-point inversion recovery and single-volume 
localization by a stimulated echo at logarithmically 
increasing TE = 20, 30, 45, 65, 95, 135 ms as described 
in (23). The agar reduces bulk motion and the T 2 /T 1 ra- 
tio, thus reducing the effect of residual transverse coher- 
ences of echo signals when compared with brain tissue. 
The agar phantom was covered by a 128 x 104 x 104 
matrix with 1 mm isotropic resolution. The nominal flip 
angle of the vendor's product FLASH sequence ("gre") 
was varied between 1° and 34°. The flip angle range was 
adjusted according to each of the four different values of 
TR (6, 12, 24, and 48 ms) to avoid crowding (6 ms: 1°- 
13°; 12 ms: 1°-18°; 24 ms: l°-24°; 48 ms: 2°-34°). Use of 
a nonselective rectangular RF pulse and a read-out band- 
width (BW) of 270 Hz/pixel resulted in a TE of 2.65 ms. 

Experiment 2: Slice-Selective Excitation 

To study the effect of slice-selective excitation in the 
vendor's FLASH sequence, bandwidth, TE, and TR were 
increased to 435 Hz, 3.2 ms, and 7 ms, respectively, 
when compared with Exp. 1. Central slices of 5 mm 
thickness were imaged on the agar phantom using the 
excitation pulse shapes ("normal" and "fast") as imple- 
mented in the vendor's "gre" FLASH sequence and com- 
pared to 3D encoding of 52 partitions across a slab of 
104 mm, that is, 2 mm partition thickness. 

Experiment 3: RF Spoiling at Varying Phase Increment 

The effect of RF spoiling was studied on a 0.1 mM solu- 
tion of MnCl 2 to smooth out the rapid variation in the de- 
pendence of T a on the phase increment (15). T a = 0.54 ± 
0.37 s and T 2 = 39.7 ± 0.2 ms were determined on two 
different occasions by single-slice 6-point inversion recov- 
ery and single-volume localization by a stimulated echo at 
logarithmically increasing TE = 20, 30, 45, 65, 95, 135, 
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200, 300 ms. In contrast to Exps. 1 and 2, a customized 3D 
FLASH sequence implemented in IDEA VB15, Siemens 
Healthcare, (21,24) was used for acquisition at TR = 7 ms, 
TE = 2.6 ms, bandwidth = 305 Hz/pixel. Based on the 
extensive simulations in (14), we chose to study two inter- 
vals of different qualitative behavior: First, the phase incre- 
ment was varied from 117.0° to 119.8° in steps of 0.4° to 
sample the rapid transition from negative to large positive 
residuals (i.e., deviations from the Ernst equation assuming 
perfect spoiling). Then, the less variable range between 47° 
and 54° phase increment was sampled in steps of 1°. This 
experiment comprises the 50° default increment of the 
"gre" sequence. For each phase increment, measurements 
were performed at VFA of 2°, 3°, 4°, 6°, 8°, 10°, 12°, 15°, 
18°, 21°, 24°, 30°, and 36°. 

Experiment 4: Human Brain In Vivo 

In vivo, an isotropic resolution of 1.25 mm and 6/8 par- 
tial Fourier acquisition was chosen to reduce the mea- 
surement time. 128 sagittal partitions of 192 x 174 pixels 
were acquired at a constant TR of 11 ms in 3:10 minutes 
per volume. The TE of 4.92 ms at a bandwidth of 200 
Hz/pixel yielded fat and water signals in phase. The 
phase increment was 50°. These settings have been used 
in a scheme to correct magnetization transfer images for 
the influence of T a and RF inhomogeneities (25). The 
flip angle was varied through a total of 16 values from 2° 
to 12° to 24° to 40° to 60° in steps of 2°, 3°, 4°, and 10°, 
respectively. A TV-weighted MP-RAGE scan (magnetiza- 
tion-prepared rapid acquisition of gradient echoes; TI = 
900 ms, a = 9°, TE = 3.2 ms, TR = 2250 ms) of 1 mm 
isotropic resolution was acquired as an individual ana- 
tomical reference. Rapid flip angle mapping was per- 
formed as described previously (26). 

Data Processing 

Mean signal and the standard deviation were assessed 
over a 10 mm region-of-interest in the center of the phan- 
toms. Kaleidagraph 3.6 for Windows (Synergy Software, 
Reading, MA) was used for spreadsheet calculations, dis- 
play, and least-squares fitting with optional weighting of 
the squared residuals. For nonlinear fitting, a modified 
Levenberg-Marquardt algorithm with numerical estima- 
tion of partial derivatives yielded the Cramer-Rao lower 
bound estimates for the standard errors. 

In vivo data were processed using the FMRIB software 
library (FSL 4.1, Center for Functional Magnetic Reso- 
nance Imaging of the Brain, University of Oxford, UK 
www.fmrib.ox.ac.uk/fsl). The MP-RAGE volume was 
aligned to the 1 mm brain template of the Montreal Neu- 
rological Institute by a rigid-body transform. The VFA 
data were then aligned to this individual reference. 3D 
region-of-interests were placed in the head of the left cau- 
date nucleus, the splenium of the corpus callosum, and 
the left lateral ventricle to evaluate signals of gray matter 
(GM), white matter (WM), and cerebro-spinal fluid (CSF), 
respectively. Maps of T t were obtained by linear regres- 
sion of signals that were unaffected by signal bias, i.e., 
that did not deviate from the straight line (JV= 8, 2° < a < 
18°). Equations 12 and 14 were calculated using custom- 



ized FSL scripts. The influence of flip angle bias was cor- 
rected post hoc on the estimated maps of T lapp (17). 

RESULTS 

Experiment 1 : Nonselective Excitation 

The phantom data followed the linear relation between 
the scaled signals well (Fig. 2a). The steeper slopes were 
observed for shorter TR, while the amplitudes were con- 
sistent within fitting errors. The standard deviations 
across the region-of-interest were taken as surrogate 
measures of image noise. Only positive error bars are 
shown to illustrate the correlation of errors along ab- 
scissa and ordinate. The T-i estimates derived from p-i 
were between 2.40 s and 2.44 s. They were consistent 
within their standard errors with the value determined 
by inversion recovery. The standard errors of the fitted 
T a decreased from 34 ms to 14 ms when increasing TR 
from 6 ms to 48 ms. This is in line with the TR-depend- 
ent decay of partially refocused echo signals. Since the 
TR's were much shorter than TV of the agar phantom, the 
abscissa was scaled by 2 TR to fit T a directly from all 
data points (Fig. 2b). This yielded a consistent value of 
2.41 s with a small standard error of 0.01 s. 

In contrast to Fig. 2a, the fitted lines at different TR 
could hardly be discerned when using the common lin- 
earization (Fig. 2c), because the slopes were close to one 
at short TR. The error bars were highly correlated. In 
particular for small flip angles (upper right), the error 
bars were almost parallel to the regression line, so the 
linear fit yielded a small error (about 10~ 4 ) for the slope, 
exp(-TR/T a ). This yielded consistent T x estimates, albeit 
with artificially low standard errors. In contrast to above, 
the standard errors increased with TR from 0.0002 ms to 
0.0067 ms. The origin region of Fig. 2c is zoomed in Fig. 
2D. The variation of the y-intercepts, A(l-exp(-TR/r a ), 
reflects the influence of TR. Again, the errors of the fit- 
ted A were smaller than in Fig. 2a by more than a factor 
of 1000. Similar errors in T a and A were observed for 
both methods after application of weights to account for 
the covariance of x- and y-errors. 

Experiment 2: Slice-Selective Excitation 

Slice selective VFA results in variable degrees of partial 
saturation across the slice profile. The signal maximum 
was shifted to higher flip angles when compared to 3D 
phase encoding (Fig. 3a). Positive signal bias was observed 
at higher flip angles since the margin of the profile yields 
a higher signal than the center. While the signal depend- 
ence on t still remotely resembled the Ernst curve, the lin- 
ear display revealed strong deviations from the 3D refer- 
ence measurement (Fig. 3b). Although the "fast" RF pulse 
(dotted) yielded about 17% higher signal, its flip angle de- 
pendence did hardly differ from that of the "normal" RF 
pulse (dashed). With 3D encoding, the slice profile added 
to the flip angle inhomogeneity and decreased T lapp in the 
off-center partitions (not shown). 

Experiment 3: RF Spoiling at Varying Increment 

Figure 4a shows the VFA measurements in MnCl 2 solu- 
tion at different RF spoiling increments. The 118.2° 
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FIG. 2. VFA measurements at different TR in 6% agar, a: Signals obtained at different TR are separated by the linear parameterization 
based on the half-angle tangent. Only positive error bars are plotted to emphasize the correlation between x errors and y errors. Close 
to the Ernst angle the resulting error is orthogonal to the line (arrow), b: VFA measurements at different TR combined into one regression 
over St/2TR. c: The traditional linear parameterization arranges the data close to the identity. Only the data obtained at TR = 6 ms and 
48 ms are shown for clarity. Data points obtained at smaller flip angles have highly correlated errors, d: Enlarged origin region of C to 
illustrate the influence of the relaxation term on the small y-intercepts. 



increment (black diamonds) conformed to a straight line 
(solid) as predicted by simulation (14). Already at 10°, 
some signals started to deviate from the linear relation- 
ship, indicating the influence of residual transverse mag- 
netization due to partial refocusing of echo pathways. 
Note that the sign of residuals was not consistent for the 
119.8° phase increment, leading to a widely deviating 
regression line (dotted). All VFA measurements at incre- 
ments between 47°, and 54° followed a similar pattern 
(Fig. 4b). Here, linearity was well obeyed up to flip 
angles of 15°, followed by increasing deviations toward 
smaller signals. 

Experiment 4: Human Brain In Vivo 

In vivo, the VFA measurement at TR =11 ms and 50° 
phase increment conformed to a straight line for nominal 
flip angles up to 15° (Fig. 5a). Data points obtained at 
higher flip angle showed the negative residues typical 



for the 50° increment and were omitted from the linear 
fit. For such flip angles up to 15°, t is practically identi- 
cal with the flip angle (16) and the correction of T a maps 
for the actual flip angle was performed post hoc as 
described in (17). RF-corrected T a values were 0.85 ± 
0.03 s for WM in the splenium, 1.29 ± 0.04 s for GM in 
the caudate head, and 4.58 ± 0.13 sec for cerebro-spinal 
fluid in the lateral ventricle. They are displayed as a 
pseudo-color overlay in Fig. 5b. The T 1 histogram 
showed the expected separation of WM and GM modes 
(Fig. 5c). 

DISCUSSION 

We derived a novel linear plot for the VFA experiment 
that has several advantages over the conventional 
method (5) that has been widely used for T x quantifica- 
tion. It provides a more intuitive display of the VFA data 
where signal amplitude and relaxation term are 
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FIG. 3. Effects of slice profile in 2D FLASH, a: Signal dependence of 2D FLASH on t with fitted Ernst curves for "fast" (dashed) and 
"normal* (dotted) selective excitation. The signal obtained with 3D encoding (solid, "normal") is shown for comparison to illustrate the 
shift of the signal maximum towards higher angles and positive residues at higher flip angle, b: Linear display of the data enhances the 
nonlinear VFA behavior of the 2D experiments, c: Conventional linear display of the 2D (with error bars) and 3D data (diamonds) with 
"normal" selective excitation. The large bias of the 2D signal is obscured, resulting in an unreasonably high linear correlation (Pearsons's 
R = 0.99999). 



represented separately as intercept and slope of the regres- 
sion line. The line segment is centered on the Ernst angle, 
which thus can easily be identified. Since the parameter- 
ization is based on the half-angle tangent, the flip angle 
range up to 180° could be exploited, in contrast to just 90° 
with the conventional method. For flip angles exceeding 
90°, however, the increasing nonlinearity of the half-angle 
tangent would require a correspondingly precise mapping 
of the flip angles. At low flip angle and TR <C T 1; how- 
ever, t and p 1 may be replaced by a and J? a TR in good 
approximation allowing for post hoc flip angle correction 

(16) . Under these conditions pertaining to FLASH imaging 

(17) , the influence of TR on the slope can be removed by 
additional scaling of the abscissa (18). Signals obtained at 
different TR are combined into a single linear regression 
that directly yields T a . This offers additional experimental 
flexibility to incorporate images measured at different TR, 



e.g., in the context of magnetization transfer experiments 
(24,25,27). 

Being a generalization of dual angle experiments 
(16,17), the linear plot yields estimates of an apparent T a 
that increases with the square of the multiplicative trans- 
mit bias. Bias correction can be performed post hoc 
using independently acquired flip angle measurements 
(26-28) or by applying a unified segmentation algorithm 
(20) to the uncorrected maps of apparent T t (21). Flip 
angle mapping constitutes an independent source of 
error that cannot be visually assessed by our suggested 
plot. The T t values determined in the splenium and cau- 
date of our single healthy volunteer, were above the av- 
erage of a cohort study applying inversion-recovery for 
T\ quantification at 3T (T a = 748 ± 64 ms in splenium; 
Tj = 1258 ± 55 ms in caudate; 29), but still fell into 
95% confidence interval. 




FIG. 4. Effects of phase increment in MnCI 2 solution, a: 117.0°-119.8° Deviations from the predicted line appear above 10°. b: 47°-54° 
Deviations from linear behavior (fitted to flip angles 2° to 8° at 52° phase increment) appear above 15°. Note, that signals obtained at 
identical flip angle fall onto lines through the origin. 
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FIG. 5. VFA measurements in vivo, a: Linear fit of signals from caudate nucleus (GM, squares, green), splenium (WM, diamonds; red), 
and lateral ventricle (cerebro-spinal fluid, diamonds, blue). Local flip angle inhomogeneities were corrected before calculating t. Devia- 
tions from the signal equation were observed for nominal flip angles above 18° and excluded, b: Pseudo-color overlay of flip angle cor- 
rected Ti values to enhance the inhomogeneities within WM and deep GM. The distribution across the cortex is degraded by partial 
volume effects at a resolution of 1.25 mm. c: Whole-brain histogram of 7"i values with color-scale of the overlay. 



The greatest practical advantage of the suggested linear 
plot is the conspicuous identification of signal bias, i.e., 
deviation from the Ernst equation representing perfect 
spoiling. This was illustrated for the partial refocusing of 
echo pathways and for slice profile effects in 2D imaging 
(22). Signal bias is difficult to identify in the conven- 
tional plot, because this imposes a large covariance of 
errors, so that even strongly deviating data points are 
projected closer to the ideal line. This resulted in unrea- 
sonably high linear correlations which may obscure infe- 
rior data quality. Such systematic errors can be even dif- 
ficult to find when fitting the nonlinear flip angle 
dependence, since the arbitrary signal amplitude is 
adjusted simultaneously. An immediate application of 
our method is to control the spoiling of RF quality of the 
FLASH sequence as implemented on one's clinical or ex- 
perimental MR system as demonstrated in this study. 

While the conventional linear plot may cover up signal 
deviation, the suggested method is highly sensitive to it 
and may thus be more prone to errors. In particular, 
incompletely spoiled echo signal became influential at flip 
angles lower than the theoretically derived 7Vw optimum 
(at 2.41 t e ). If suitable correction techniques (14) are to be 
applied, one still can pursue SNR optimization. However, 
if such processing is not available, it is recommended to 
sacrifice SNR in favour of reducing systematic errors in 
the fitted parameter maps (as done in Exp. 4). 

The contributions of spin echoes and stimulated echoes 
vanish in the limit of a — > 0, because they depend on 
higher order powers of the flip angle. Since the signal bias 
is negative at most phase increments, it will be enhanced 
when approximating t by a. Spoiling gradient areas of 
280-450 mT ms / m as recommended by (15) to reach the 
diffusion limit can only be achieved on clinical MRI sys- 
tems by increasing the spoiler duration and thus TR. Sig- 
nal offsets due to Rician noise at low SNR and motion (22) 
will cause additional deviations. The negative residues 
observed at very low flip angle (l°-2°) we tentatively as- 
cribe to additional saturation due to noise leaking into the 
sample from the unblanked RF power amplifier. Measure- 



ments of rapidly relaxing water moieties at ultra-short TE 
will also deviate from the ideal flip angle dependence (30). 

A full mathematical account of the method was 
derived, as a generalization of the dual angle approach 
(16). The signal equation of partial saturation was trans- 
formed into an algebraic function of t obtained by the 
half-angle tangent substitution of trigonometric func- 
tions. The relaxation term was cast into a single expres- 
sion (px) by a similar transformation, yielding a low- 
order rational function in t and p a . This facilitated the 
use of computer algebra for deriving the regression for- 
mulae (Eqs. 12 and 14) and the propagation of image 
noise to the TV and amplitude maps. As an important 
result we found that there was no common optimal set- 
ting for maximal SNR in both TV and amplitude maps. 
Theoretical noise optimization for N measurements 
demanded a clustering of flip angles; particularly a repe- 
tition of the dual angle scheme to estimate T a from an 
even N. This is in contradiction to results obtained with 
a genetic algorithm for the conventional linear parame- 
terization (8). This may be due to behavior of the genetic 
algorithm or to differences in the scaling of signal noise 
which may affect linear regression analysis (10). Suitable 
weights for the residues have been derived for the sug- 
gested methods in the theory section, but as a rule these 
require an a priori estimate of T a . However, the effects of 
signal bias most often exceeded those of signal noise. 

The immediate application of our method refers to the 
flip angle dependence of the residual signal deviations. 
The main rationale for acquiring signals yielding equi- 
distantly spaced points on the linear plot — even at the 
cost of reduced SNR — is to define the straight line, 
against which signal bias at low and high flip angles can 
be identified and excluded from an iteration of the 
regression. This applies also to patient motion, which 
can corrupt singular volumes of the VFA series. The in 
vivo experiment was performed with a phase increment 
of 50° as implemented in the MR scanner's FLASH 
sequences and previously recommended for dual angle 
T-i mapping (14). The flip angle dependence of the signal 
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conformed to the Ernst equation for rather small flip 
angles (up to 15°-20°). These may be well lower than the 
higher flip angle in a noise-optimized dual angle proto- 
col. The flip angles for noise optimized T a (Table 1) cor- 
respond to 85.3% and 14.7% of the amplitude level in 
S/t. The latter signal may be well affected by bias, which 
can be corrected by suitable calibrated models (14). 
Alternatively, deviating data points at high flip angle 
data have to be identified and subsequently excluded 
from the analysis. The theoretical optimization does not 
take into account signal bias, but still confirms the intui- 
tive insight that the precision of the amplitude benefits 
from acquisitions at low flip angles. High resolution am- 
plitude maps have been shown to be useful for depicting 
deep brain nuclei (31). In fact, amplitude mapping corre- 
sponds to correcting the measured signals for varying 
relaxation effects; which is achieved here by extrapola- 
tion of S(t)/t to t = 0. The conventional method, on the 
other hand, places data points closer to the y-axis as the 
flip angle approaches 90°. The contributions of partially 
refocused echo pathways will thus influence the estima- 
tion of y-intercept, and hence the amplitude maps. 

Our study imposes practical consequences for the 
design of a VFA protocol. Upon choosing the TR and the 
number of measurements, data points should be evenly 
distributed along the regression line to fully exploit the 
potential of the suggested linear plot to identify bias. To 
account for local flip angle bias, the corresponding data 
points are centered on the apparent Ernst angle, that is, 
prior to flip angle correction. This condition cannot be 
met across the whole brain. The median of the apparent t e 
can be a suitable compromise to determine the array of flip 
angles from Eq. 11. This scheme may require values of 
fractional degrees. If T a is to be estimated without correc- 
tion other than the indispensable RF correction, it should 
be avoided to acquire biased signals at very low flip angles 
(=2°) and those producing partially refocused transverse 
coherences at higher flip angles. An iterative exclusion of 
such points becomes possible with our scheme. 
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